Magnetic properties of Mn-doped Ge46 and Ba8Ge46 clathrates 
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We present a detailed study of the magnetic properties of unique cluster assembled solids namely 
Mn doped Ge46 and BagGe46 clathrates using density functional theory. We find that ferromagnetic 
(FM) ground states may be realized in both the compounds when doped with Mn. In Mn2Ge44, 
ferromagnetism is driven by hybridization induced negative exchange splitting, a generic mechanism 
operating in many diluted magnetic semiconductors. However, for Mn-doped BagGe46 clathrates 
incorporation of conduction electrons via Ba encapsulation results in RKKY-like magnetic interac- 
tions between the Mn ions. We show that our results are consistent with the major experimental 
observations for this system. 



Cluster assembled solids formed by large cages of Si 
and Ge, the so called clathrates have drawn consider- 
able attention in the recent years because of the wealth 
of physical properties displayed by them. These periodic 
arrangement of nano-sized cages may be semiconduct- 
ing with energy gap larger in comparison to that of the 
standard diamond form of Si and Gc.^ Alkali/alkaline 
earth encapsulated Si and Ge clathrates are metals.^ 
Some of them are found to be superconducting'^ and are 
also suggested to be excellent candidate for thermoelec- 
tric applications.'* The germanium clathrates were also 
found to be ferromagnetic upon doping with transition 
metals (TMs).^^^ Ferromagnetism was also reported for 
Eu4Sr4GaigGe3o^ and EugGaieGeao^ clathrates. The 
discovery of magnetism in Ge-based clathrates have 
added a novel functionality to these cluster-assembled 
solids that may find possible application in magnetic de- 
vices. 

Ge46 is a type I clathrate with simple cubic lattice, 
where the Ge atoms form closed cage-like structures of 
Ge2o and Ge24. The unit cell consists of two 20 atom 
(Ge2o) cages and six 24 atom (Ge24) cages. Among the 
two Ge2o cages one is placed at the corner of the cube, 
while the other, rotated 90° with respect to the first is 
placed at the body center of the cube. Each Ge in Ge2o 
cage is bonded to three other atoms within the same cage. 
The four fold coordination is completed in two steps. The 
eight Ge atoms in each Ge2o cage form bonds along the 
eight (111) directions with the other Ge2o unit and the 
co-ordination of the remaining 12 atoms in each Ge2o 
cage is taken care by adding six Ge atoms at the inter- 
stitial 6c sites in the primitive unit cell. The resulting 
primitive unit cell has 46 Ge atoms (space group Pm3n) 
where each Ge is tetrahedrally bonded and therefore the 
material is expected to be semiconducting. Although 
semiconducting pure Ge clathrates can be experimentally 
synthesized,^" they however prefer to encapsulate alkali 
and alkaline earth metals like Na, K, and Ba, which oc- 
cupy the centers of the Ge2o and Ge24 cages (see Fig. 1). 
Such encapsulation induces a transition from semicon- 
ducting phase to a metallic phase. 

The interstitial 6c sites (highlighted in Fig. 1) are of 
particular importance as the doping of Mn atoms ran- 




FIG. 1: (Color online) Structure of BagGe46 clathrate show- 
ing the simple cubic lattice with Ge2o and Ge24 cages. Ge 
atoms at the crystallographic 6c sites are highlighted. All the 
Ge atoms in the unit cell are not shown for clarity. 



domly at these 6c sites of Ba8Ge46 clathrate was reported 
to yield ferromagnetism with saturation magnetic mo- 
ment per Mn atom to be ^0.8 iib and the ferromagnetic 
Curie temperature to be ^10 K.^ The distance between 
the dopant Mn atoms at the 6c sites was also found to 
be quite large, ^ so it was speculated that the interaction 
between the dopants may be Ruderman-Kittel-Kasuya- 
Yosida (RKKY)-like. Following this report, Yang et al.^'^ 
have investigated this system in the framework of ah ini- 
tio density functional calculations. For a pair of Mn 
atoms to be substituted at the 6c sites in the unit cell, if 
the position of one Mn is fixed then the other can occupy 
any of the remaining five sites. Out of these five neigh- 
bors of Mn, one is at a distance a/2 (configuration I) while 
the other four are at a distance •\/6a/4(configuration II), 
where a is the lattice constant. The calculation by Yang 
et al. in the framework of LDA yielded the value of 



the magnetic moment per Mn atom to be 0.77 /is and 
0.42 hb when the two Mn atoms are doped in configura- 
tion I and configuration II respectively^^. Although the 
former value of the magnetic moment agrees well with 
the experiment but the latter is found to be substan- 
tially less. Further, Yang et al.^'' did not explore the 
possibility of antiferromagnetic (AFM) coupling between 
a pair of Mn atoms necessary for RKKY like interaction. 
In view of the above, in this letter we have studied the 
magnetic properties of Mn doped Ba8Ge46 and pristine 
Ge4g clathrates in some details using ah initio density 
functional calculations. Our calculations reveal that fer- 
romagnetism may be realized in both the systems but the 
origin of ferromagnetism in guest-free Mn2Ge44 system is 
markedly different from metal encapsulated BagMn2Ge44 
clathrates. 

The electronic structure and total energy calculations 
presented in this paper are performed using ab initio 
density functional theory (DFT) as implemented in the 
Vienna Ab-initio Simulation Package (VASP).^'''^^'' The 
electron-ion interaction in the core and valence part are 
treated within the projector augmented wave (PAW) 
method^^ along with plane wave basis set. We have em- 
ployed generalized gradient approximation (GGA) due to 
Perdue-Burke-Ernzerhof (PBE)^'^ to treat the exchange 
and correlation. The localized Mn d-states are treated in 
the framework of GGA-I-U method, where calculations 
are done for several values of U in the range 2 — 7 eV 
and J=l eV. Atomic positions were relaxed to minimize 
the Hellman-Feynman forces on each atom with the tol- 
erance value of 10~^ eV/A. The optimum values of the 
energy cutoff and the size of the fc-point mesh is found to 
be 650 eV and 4x4x4 respectively and was accordingly 
employed in our calculation. 

To begin with we have computed the electronic struc- 
ture of pristine Ge46 and Ba8Ge46 clathrates without 
doping. Our optimized lattice constants for Ge46 and 
Ba8Ge46 clathrates are 10.78 A and 11.01 A respec- 
tively. The density of states(DOS) corresponding to these 
clathrates are displayed in Fig. 2. We see from Fig. 2(a) 
that Ge46 is a semiconductor with a band gap of 1.19 eV. 
The valence band DOS of the pristine Ge46 (Fig. 2(a), 
from left to right) exhibits three major parts as reported 
earlier^° that may be assigned to a s-like region, a sp- 
hybridized region, and a p-like region, with a character- 
istic s-p gap. The origin of this s-p gap has been at- 
tributed to the five ring patterns (pentagons) of the Ge 
atoms. We have displayed the total DOS as well as 
the partial density of states (pDOS) for Ba encapsulated 
Ge46 clathrate in Fig. 2(b) and Fig. 2(c), (d), and (e) 
respectively. Upon Ba encapsulation the additional 16 
valence electrons are accommodated in the conduction 
band and the system becomes metallic. While the DOS 
of the valence band of Ba8Ge46 is very similar to pris- 
tine Ge46, the conduction band DOS shows modification 
upon the inclusion of metal atoms. The pDOS plot of 
Ba8Ge46 clathrates (see Fig. 2(c), 2(d), and 2(e)) reveal 
that the conduction band has considerable admixture of 
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FIG. 2: (Color online) The density of states corresponding 
to Ge46 and BagGe46 clathrates. Subfigure (a) and (b) show 
the total density of states for Ge46 and BagGe46 clathrates 
respectively. Ge-s, Ge-p and Ba-s -I- Ba-d states correspond- 
ing to the Ba8Ge46 clathrate are depicted in subfigures (c), 
(d) and (e) respectively. 



Ge-s and p states with Ba-s and Ba-d states with sub- 
stantial contribution from the latter. These electrons are 
expected to be dclocalized in the entire system and act 
as a conduction electron cloud. 

Next we have considered doping of a pair of Mn 
atoms at the 6c sites of the semiconducting pristine Ge46 
clathrate. In order to understand the magnetic proper- 
ties of the doped system we have calculated the magnetic 
exchange interactions as the energy difference between 
antiparallel (AFM) and parallel (FM) arrangements of 
the magnetic moments of a pair of Mn ions for config- 
urations I and II for several values of U ranging from 
U=0 (GGA) to U=7 eV. Our calculations reveal that the 
exchange interaction is always ferromagnetic in both the 
configurations for all values of U except for U=7.0 eV (see 
table SI-1).^^ These results suggest that Mn doped Ge46 
may be ferromagnetic. For further insights on the ori- 
gin of magnetic ordering in Mn2Ge44 we have displayed 
in Fig. 3 the spin polarized density of states (DOS) for 
Mn2Ge44 with parallel alignment of the Mn spins (FM) 
in configurations I for the representative case, U=4 eV 
and J=l eV. The DOS indicates that the system is a 
half- metallic ferromagnet, analogous to diluted magnetic 
semiconducting (DMS) materials like Mn-doped GaAs. 
When a Ge is substituted by a Mn having seven valence 
electrons, four electrons are utilized to saturate the Ge 
dangling bonds and the remaining three electrons are re- 
sponsible for the magnetism. A comparison of the to- 
tal DOS (Fig. 3(a)) and the pDOS for Mn-d (Fig. 3(b)) 
and Ge-p (Fig. 3(c)) suggests that the Mn-d states are 
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FIG. 3: (Color online) The spin polarized GGA-hU (U=4 
eV and J=l eV) density of states for Mn2Ge44, with paral- 
lel alignment of Mn spins (FM) in configuration I. Subfigure 
(a) shows the total DOS, while subfigures (b) and (c) indi- 
cate the projected DOS corresponding to Mn-3d and Ge-4p 
states. Subfigure (d) represents the partial density of states 
corresponding to one dopant Mn atom and its surrounding Ge 
atoms, whereas the inset shows 2-dimensional charge density 
for the up spin valence band hole states. 



not only strongly exchange split but also Mn-d spin up 
states are fully occupied and the Mn-d spin down states 
are completely empty confirming that Mn is in high spin 
d^ configuration. So the introduction of Mn in pristine 
Ge46 gives rise to an acceptor in the Ge -p manifold and 
the hole produced by the acceptor (see Fig. 3(b)) is ex- 
pected to interact with the localized Mn-d states to me- 
diate fcrromagnetism. We find from Fig. 3(b) that the 
exchange split Mn 3d states in the ferromagnetic con- 
figuration strongly hybridizes with the Ge-p states. As a 
result, the spin up Ge-p bands shift to higher energy while 
the spin down Ge-p bands make an opposite shift as indi- 
cated by the arrows in Fig. 3(b) and therefore an energy 
gain is obtained by transferring electrons from the Ge 
spin up states to the Ge spin down ones. The exchange 
splitting of the p-states induced by this hybridization are 
opposite to that of the Mn-d states (i.e. hybridization in- 
duced negative exchange splitting) resulting in the calcu- 
lated magnetic moment of Mn and Ge oppositely aligned 
and the total magnetic moment to be 6 fiB per unit 
cell. This generic mechanism also responsible for fcr- 
romagnetism in Sr2FeMo06^'^"^^ and Mn doped GaAs^^ 
is further corroborated by the plot of pDOS (Fig 3(d)) 
and charge density of the valence band hole states in the 
spin up channel, shown in Fig. 3(d) (inset). From these 
plots we gather that there is strong hybridization of the 



Mn— d states with the Ge—p states, thereby accounting 
for fcrromagnetism in Mn doped Ge46 . In the presence of 
strong hybridization between the Mn-d and Ge-p states 
this novel mechanism of fcrromagnetism underscores the 
competing antiferromagnetic ordering of Mn spins via su- 
per exchange. 

Finally we have considered the doping of Mn in 
Ba8Ge46. We have first assumed substitutional doping of 
two Mn atoms into the Ba8Ge46 clathrate at the 6c sites. 
The equilibrium lattice constant and the atomic posi- 
tions are calculated for the Mn doped at the 6c sites both 
in configuration I and configuration II for the ferromag- 
netic and antiferromagnetic alignment of Mn spins. The 
equilibrium lattice constant is calculated to be 11.08 A, 
which is an overestimation from the experimentally re- 
ported value by 3.6 % possibly due to the use of GGA 
for the exchange-correlation. •^'^ We have also calculated 
the exchange interaction and the magnetic moment per 
unit cell (in the FM state) in configuration I and con- 
figuration II for values of U ranging from U=0 (GGA) 
to U=7 eV (see Table SI-2).22 We find that in configura- 
tion I the exchange interaction is antiferromagnetic while 
in the configuration II, it is ferromagnetic, independent 
of the chosen values of U (see Table SI-2).^^ This is in 
sharp contrast to the case of Mn doped Ge46. The change 
of the magnitude as well as the sign of the exchange in- 
teraction with distance points to the fact that the RKKY 
model may be applicable here. It is well known that for a 
fixed concentration of the magnetic ions, the RKKY in- 
teraction depends only on the distance between the mag- 
netic ions in the presence of the delocalized conduction 
electron cloud. So in order to extract the characteristic 
distance dependence of the RKKY interaction, we have 
also considered doping configurations, where one Mn is 
at the 6c site but the second one is in other possible 
framework sites as the conduction electron cloud spread- 
ing the entire clathrate network mediating the exchange 
interaction will possibly not distinguish among the dif- 
ferent framework sites. The results of our calculations 
for the exchange interaction for the representative case 
U=4.0eV and J=1.0eV are summarized in Fig. 4. We 
have fitted the calculated exchange interactions with the 
RKKY expression for the exchange interaction strength 
J(r) given by 

J(r) = Const. X [sin(2fcFr) - 2kFr C08{2kFr)]/r^ , (1) 

where r is the distance between dopant magnetic atoms 
and kp is the Fermi wave vector corresponding to the av- 
erage electron density.^^ Fig. 4 shows excellent agreement 
of the fit with RKKY model with the fitting parameter 
kp = 1.05 A^^, confirming RKKY mechanism to be op- 
erative in Mn doped Ba8Ge46 clathrates. Further the plot 
of the charge density shown in Fig. 4(inset) reveal that 
it is delocalized in the entire system and the exchange 
interaction is mediated by the conduction electron cloud. 

We shall now compare our calculated results with the 
available experimental data. In view of the large average 
distance (d) between Mn atoms located at the 6c sites 
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FIG. 4: (Color online) The variation of magnetic exchange 
interaction strength with distance between the Mn atoms has 
been plotted and fitted to the RKKY expression with kp = 
f.05 A~^. Inset shows the 2-dimensional charge density plot 
on a plane containing one dopant Mn atom. 

(d=8.5A for lattice constant a = 10.689 A ) reported 
for this system^ we have extrapolated the RKKY plot to 
obtain the exchange interaction between the Mn ions lo- 
cated at the 6c sites other than configuration I and config- 
uration II. Fig. 4 reveal that the exchange interaction is 
not only ferromagnetic for configuration II (d= 6.79 A for 
theoretically calculated lattice constant a=11.08 A) but 
also when the separation between a pair of Mn atoms at 
the 6c site is d=10.36 A. The large interval between the 
Mn atoms reported in the experimental work^ possibly 
refers to the latter two configurations where the exchange 
interaction is always ferromagnetic. The small saturation 



moment (0.8 ^b) obtained experimentally finds a natural 
explanation in such an RKKY scenario. The Mn loca- 
tions are random due to competitive total energies (either 
FM or AFM) and the Mn atoms that prefer the AFM 
state do not contribute to the magnetic moment and the 
moment is due to a small fraction of Mn atoms that are 
magnetically active in the FM state. The weak exchange 
interaction^^ for Mn doped Ba8Ge46, also suggests a low 
Curie temperature in BagMn2Ge44, in agreement with 
the experiment. 

In conclusion, we have studied the electronic structure 
of Ge4g and Ba8Ge4e clathrates doped with Mn from 
ab initio density functional calculations. We find that 
ferromagnetic ground state may be realized in both the 
compounds. The origin of ferromagnetism in Mn2Ge44 is 
driven by hybridization induced negative exchange split- 
ting while it is RKKY-like for Mn doped Ba8Ge46. The 
origin of the two different mechanism may be traced 
back to the electronic structure of these systems and 
limit of validity of the RKKY model. The RKKY limit, 
<< 1, where is the exchange splitting of the 
host band and Ep is the Fermi energy is not satisfied 
for the half-metallic system Mn2Ge44 due to complete 
spin polarization, resulting in E^ > Ep.^^ However in- 
corporation of the conduction electrons in Ba8Mn2Ge44 
upon Ba encapsulation makes the system metallic and 
protects the RKKY limit. The RKKY-like scenario pre- 
dicted for Ba8Mn2Ge44 is also consistent with the major 
experimental observations for this system. 
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Supplementciry Information 

TABLE SI-1: Energy difference between antiferromagnetic and ferromagnetic states and the magnetic moments per unit cell 
(in the ferromagnetic state) for Mn2Ge44 in both configurations and different values of Hubbard U parameter. 



U (eV) 


Configuration I 

Eafm — EpM magnetic moment 
(meV) (/is per cell) 


Configuration 11 

Eafm — Efm magnetic moment 
(meV) {fiB per cell) 


0.0 


186.3 


6.0 


151.1 


6.0 


2.0 


182.5 


6.0 


145.2 


6.0 


4.0 


171.6 


6.0 


134.3 


6.0 


5.0 


21.5 


7.6 


14.4 


7.7 


6.0 


10.3 


8.2 


5.2 


8.2 


7.0 


-6.1 


8.5 


-8.5 


8.7 



TABLE SI-2: Energy difference between antiferromagnetic and ferromagnetic states and the magnetic moments per unit cell 
(in the ferromagnetic state) for Ba8Mn2Ge44 in both configurations and different values of Hubbard U parameter. 



U (eV) 


Configuration I 
Eafm — Efm magnetic moment 
(meV) {hb per cell) 


Confi^ 
Eafm — Efm 
(meV) 


juration II 
magnetic moment 
{jiB per cell) 


0.0 


-22.8 


6.1 


22.7 


6.1 


4.0 


-54.2 


6.3 


19.9 


6.4 


5.0 


-56.4 


7.5 


12.2 


7.6 


7.0 


-142.0 


9.1 


9.3 


9.2 



